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This is a summarising investigation of the events of the phase transition of the first order that occur in the 
critical region below the liquid-gas critical point. The grand partition function has been completely integrated 
in the phase-space of the collective variables. The basic density measure is the quartic one. It has the form of 
the exponent function with the first, second, third and fourth degree of the collective variables. The problem 
has been reduced to the Ising model in an external field, the role of which is played by the generalised chemical 
potential p*. The line n* (rf) = 0, where q is the density, is the line of the phase transition. We consider the 
isothermal compression of the gas till the point where the phase transition on the line p* ( r /) = 0 is reached. 
When the path of the pressing reaches the line p* = 0 in the gas medium, a droplet of liquid springs up. 
The work for its formation is obtained, the surface-tension energy is calculated. On the line p* = 0 we have a 
two-phase system: the gas and the liquid (the droplet) one. The equality of the gas and of the liquid chemical 
potentials is proved. The process of pressing is going on. But the pressure inside the system has stopped, 
two fixed densities have arisen: one for the gas-phase rjg = 77 c (1 - di2 ) and the other for the liquid-phase 
= 77c(l + dl 2) (symmetrically to the rectlinear diameter), where q c = 0.13044 is the critical density. Starting 
from that moment the external pressure works as a latent work of pressure. Its value is obtained. As a result, 
the gas-phase disappears and the whole system turns into liquid. The jump of the density is equal to q c d, 
where d = \/D/2G ~ t vI2 . D and G are coefficients of the Hamiltonian in the last cell connected with the 
renormalisation-group symmetry. The equation of state is written. 

Key words: critical point, phase transition of the first order, grand partition function, quartic density measure, 
collective variables 

PACS: 05.70.Jk, 64.70.F-, 64.60.F- 


1. Introduction 

The present research follows our previous consideration concerning the liquid-gas critical point. In 
this work we continue the investigations of the system of interacting particles in the vicinity of the critical 
point (THU. The great partition function was calculated. The potential of the interaction between particles 
consists of the sum of the van der Waals interaction and of the short-range interaction between the hard- 
spheres. We used the collective variables (CV) method. The Jacobian of the transition from the Cartesian- 
coordinates space to the space of the collective variables was averaged with the statistical weight of a 
distribution function of the hard-spheres. 

The problem was reduced to the Ising model in an external field, presented by the generalised chem¬ 
ical potential p*. 

As a basic density measure, the quartic distribution function is used. It has the form of the expo¬ 
nent function with the first, second, third and fourth degree of the collective variables. Detailed calcu¬ 
lations were performed in 0 for the case of t 3= 0, where r = (T- T C )IT C . Therein, the coordinates of 
the critical point T c and rj c were obtained, where 7 C is the critical temperature, rj c is the critical density 
q c = ( NIV) C (ttct 3 /6), a is the diameter of a hard sphere. 
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Table 1. The critical indexes for the quartic and sextic density measures in the CV method (5). 


Critical parameters 

V 

a 

P 

7 

Approximation p 4 
Approximation p 6 

0.605 

0.637 

0.18S 

0.088 

0.303 

0.319 

1.210 

1.275 


The critical indexes correspond to that of the Ising model. Their values, the sizes of the particles, the 
critical temperatures and densities are given in tables[T|and[2](T|- 

The first results of the investigations in the area of temperatures below the critical point, rsO, are 
given in GH3. It has been shown that the events connected with the phase transition of the first order 
from the gas phase to the liquid phase (and inversely) are described by the integral over the variable 
po for k = 0. This integral embraces the results of the integration over all collective variables with the 
exception of the variable po- 

The collective variable po is connected with the distribution of the number of particles. The integra¬ 
tion over the variable po was performed in the steepest descent method. It is shown that the discriminant 
Q of the equation concerning the maximum of integrand can take positive or negative values, or be equal 
to zero. When the discriminant Q is positive we have pure monophasic states of the gas or of the liquid. 
In the area Q < 0, the maximization problem possesses three real roots, 

„ av , <p + 2nn 

p™ = b cos—— -, n — 0,1,2, O^cp^n. 

However, only two of them, actually for n - 0 and n = 1, correspond to a maximum of the integrand. 
The areas of the principal and of the secondary maxima were determined. It was shown that there is a 
break between the end of the principal maximum for the gas phase and the beginning of the principal 
maximum for the liquid phase. Those were the main results of the papers (l|[2). 

In this work we would like to pass the whole route of the phase transition of the first order. The work 
consists of the introduction and seven sections. In Introduction we define the problem and give the main 
results of the mentioned works. In section [2] the description of the monophasic states is given, which 
corresponds to the positive values of the discriminant Q. Section [5] is devoted to the coordinate form 
of the reference system. We work within the framework of the quartic density measure. The reference 
system is included in our problem through the cumulants 9Jli, ')J( 2 ,11(3 and IIi 4 . In accordance with the 
accuracy of the consideration of the whole problem, it is needed to take a suitable expression for the 
pressure Po, the chemical potential po and the free energy Fq of the reference system. 

Sectionals an essential one in this work. The domain of the phase transition is considered. Here, the 
discriminant is negative Q < 0. In this region we have two maxima — the principal and the secondary one. 
The principal maxima correspond to the pure gas or to the pure liquid phases. The secondary maxima 
are connected with the stable density fluctuations, which are either of the liquid type in the gas phase, or 
of the gas type in the liquid phase. 


Table 2. The critical temperatures T c and the effective hard sphere diameter a for some systems. 


System 

P c , atm 

P c , atm 

a 

r c , °C 
(exp.) 

T c , °C 

m 

cr 0 , A 
(exp.) 

(7, A 

oIoq 

elk B , K 

CO-CO 

34.532 

37.92 

-140.23 

-138.46 

3.76 

3.37 

0.898 

100.2 

Ar-Ar 

47.964 

51.75 

-122.65 

-123.27 

3.405 

3.14 

0.922 

119.8 

Kr-Kr 

54.182 

56.76 

-63.1 

-67.84 

3.6 

3.367 

0.935 

171 

Xe-Xe 

57.537 

65.88 

16.62 

16.84 

4.1 

3.71 

0.905 

221 

O 2 -O 2 

49.77 

54.76 

-118.84 

-110.8 

3.58 

3.18 

0.89 

117.5 

N 2 -N 2 

33.55 

44.08 

-147.05 

-150.02 

3.698 

3.365 

0.91 

95.05 
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Under consideration is the equation for the envelope curve of two maxima, the principal and the 
secondary ones: 


PV , 

-q In— p 0 - 0, 


dlnHpo 

<5/lp 


= N. 


Having opened the second equation we get the probabilities of the system to be in the states correspond¬ 
ing to the principal or to the secondary maxima. We have got analytical expressions for the probabilities 
to be in the gas or in the liquid phases. 

The central considerations concern the events along the line p* = 0. Here, when the gas system is 
pressed, a droplet of liquid appears inside it, that is, a new phase arises — the phase of liquid. In the ap¬ 
pearance of the droplet lies the essence of the phase transition of the first order (as a result of isothermal 
pressure on the gas system). The work for the formation of a droplet — the surface-tension energy — was 
calculated. 

A fundamental quantity A is introduced. It describes an inverse to the usual order parameter, char¬ 
acterizing the weight of the pure monophase state in the system. On the rectangle vertex of the Q - 0, 
the weight is equal to unity for the gas phase and zero for the liquid phase or, on the contrary, dependent 
on which of them the maternal is. The quantity A is equal to zero on the line of the phase transition. 
This means that on that line the gas phase (maternal) and the liquid phase (in a droplet-form) have equal 
statistical weights. In its initial determination A was a function of the reference system cumulants 3)12. 
9K 3 ,3)1 4 and, therefore, was a function of the compressibility coefficients and its derivatives. The equality 
A = 0 means that on the line p* = 0, the compressibility coefficient equals zero. We have some metastable 
state — the state of phase transition of the first order. 

In section [6] the equality of chemical potential of the gas phase and of the liquid phase on the line 
p* = 0 is shown. 

In section [7] the gas-condensation process is analyzed. An analytical expression for the surface-ten¬ 
sion energy is obtained. The latent work Ai at in the isothermal pressure process (an analogue to the latent 
heat in an isobaric process), which acts along the line p* = 0 during the gas-liquid phase transition, was 
obtained. 

And, in the end, in section [8] the equation of state was obtained. The equation for the isotherm P = 
p (??)T=const has three parts: 


1) a monophasic gas branch: P = P(p,r) for densities p ^ p g = p c (l - d/2), d - \JD/2G (D and G are 
coefficients in Hamiltonian at the second and at the fourth degree of po); 


2) a horizontal line p* (p) = 0 for the two-phase system, here P = 0, and p g p ^ pi, pi — p c (l + d/2). 
The transition of gas into liquid is takes place at the expense of the latent work Ai at . acting from 
without; 


3) the liquid-branch for p > pi = p c (l + d/2). 


Thus, during the gas-condensation process we get a jump of the density of the amount of p! - p g = 
T] c d ~ t vI2 . As we see, the magnitudes of the densities of the beginning and of the end of the transition 
process from the gas state to the liquid state are located symmetrically to the rectilinear diameter. 

When the system approaches to the critical point, the specific surface-tension energy disappears as 
t 5/2v , because D ~ t 2v , G ~ x v , and the jump of the density p is proportional to t v/2 . 

In the preamble we stressed that within the framework of the Gibbsian statistics one can describe the 
phase transition of the first order @. 

As for other works, we benefited from the monographs of R. Balescu Q, L.D. Landau and E.M. Lifshits 
0D> J.-P. Hansen and I.R. McDonald (9j as well as from flOHT2l . 

In this research we made a number of assumptions, especially concerning the reference system, the 
dependence of the coefficients D and G on the density and others. In fact, by means of the method of 
collective variables it is possible to describe the critical point theory of the liquid-gas system with a more 
sufficient accuracy. Now we get to the theme. 
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1.1. How it should look like from the very start 

Let a big volume Vo in the form of a sphere with the radius Rq, Vo = % Rq and with the number of 
particles IV, N = no Vq be picked out in a gas medium with the density no, pressure P at the temperature T. 

We consider the phase transition of the first order as a result of the quasistatic pressure on the gas 
sphere. The process takes place at temperatures close to the critical point T c . The dimensionless temper¬ 
ature r is introduced r = (T— T C )IT C , and |r[ 0,01, r < 0. The pressure inside the sphere of radius Rq is 

equal to P + 8P, where 8P is the additional pressure, connected with the spherical form of the volume V. 
The gas energy inside the sphere: 

g 0 = (P + 8P)V 0 -aS 0 , (1.1) 

where So is the area of the sphere, So = 4 tiRq, a is the surface-tension coefficient. The minimum of the 
difference So - PVq means that 8{So - PVq)/8R - 0 or 

„ 2 a 

SPo = — ■ (1.2) 


Such a picture may be useful to remember when discussing the isothermal phase transition in the gas- 
liquid system illustrated in this work. 

We consider the grand partition function 


OO yN 

S =I^* 

N= 0 JV ' 


(1.3) 


where N is the number of particles, z N is the activity 


z 


N 


ltnk B T\ 312 1 

V 27T J h 3 


N 

expiP/dN), 


(1.4) 


m is the mass of a particle, ko is Bolzman’s constant, T is temperature, h is Plank’s constant, /3 = [k^ T)~ l , 
H is the chemical potential, Z N is the configuration integral of N particles (TJGO. The collective-variables 
method is used {p^} J3JJ5 ] to calculate E. The Jacobian of transition from the Cartesian coordinates to the 
collective variables was calculated with the statistical weight of the system of hard cores. The latter was 
taken as a reference system. As a result, for E we get the expression (l][2): 


S = HoEiH po , (1.5) 

where Ho is the partition function of the reference system (RS), 

lnE 0 = piio(N) - /3F 0 = , (1.6) 

po is the chemical potential of RS, Fq is free energy of RS. 

The quantity Hi is the partition function (T), which includes all interaction effects, with the exception 
of the effects connected with the formation of a new phase, i.e., the liquid phase within the maternal 
sphere of gas in the volume Vo'. 

Hi = Egexp [-/)(F C r +-Pigr)] ■ d-7) 

Here, E g is the partition function in which the effects appearing after integration over the collective vari¬ 
ables pk for k> B are taken into account, where B is the point of the first zero of the Fourier transform of 
the attraction potential; the integration is performed with the Gaussian density measure; Fcr + Figr are 
the free energies arising after integration over p^ with 0 < k sS: T'cr in the critical regime, where the 
renormalization group (RG) symmetry between the block-Hamiltonians exists and 2nn T ! a, h k s= n T 
is the number and a nj is the periods of the block-lattice, on which, given r, the renorm-group symmetry 
between the coefficients of the block-Hamiltonians terminates; the integration is performed in the quartic 
density measure (3j|5); Figr appears as a result of the integration over pk performed here in an inverse- 
Gaussian regime (3j for 0 < k 2jrn T l a rli . We suppose that all expressions in Hi are known, and that all 
of them are fixed quantities within the temperature region r < 0.01. In this work, we consider only one 
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part, E po , of the E. It is connected with the density fluctuations of the number of particles inside of the 
sphere of the radius Rq. The phenomena of the phase transition of the first order at t < 0 are described 
by the expression E po . The chemical potential p of the whole system is only within E po . As a result of 
calculations carried out in |T|[2), the initial form of the partition function E po for integration over po is as 
follows: 

E po = exp[p*(l-A)iv] J exp[NE(po)]dpo, ( 1 . 8 ) 

where p’ is the generalized chemical potential of the system^] 

p* =£( J u-Mo) + 'Ma(0)l(l-A), (1.9) 


po is the chemical potential of the RS, 




|9JJ 4 | ; 


A=-|^J( 2 +-TO 3 f 


( 1 . 10 ) 


( 1 . 11 ) 


SR 2 , ilK. 3 , 11(4 are the cumulants of the Jacobian of transition CHS); at small densities, A is a linear function 
of density 77 , see figure[ 2 ] 

/Ar<t>(k)) 

1 ' ( 1 . 12 ) 

(1.13) 

(1.14) 


a( 0 ) = 

[VO 

where 3>(k) is the Fourier transform of the attraction potential; 

Eipo) = p* p 0 + Dpl~Gp 4 0 , 


D - Dq\t\ 


2v 


G=Go|t|' 


are coefficients of the limiting block-structure Hamiltonian H„ t at the variable po, £>o - 1.19; Go = 1.67; 
v — Ins* /ln£i = 0.605 is the correlation length critical exponent, £1 is the greater of two eigenvalues 
of the matrix for linearized recurrent relations of Wilson’s type [ 3}|5] 11314161 for coefficients of block- 
structure Hamiltonians, s* is the optimal parameter in the scale transformations s* - 3.58.. ., (5*) < £1 < 
(s*) 1 2 3 . We work in a critical region of temperatures and densitie^] r < t* ,77 < 77 *. For the limit interval of 
the temperature t* we have the expression: 0) 


ClT* 

r* +700(0)77 X 


1 -- 


772 




1/2 


(1.15) 


Here, r* is one of the fixed-point coordinates, a 2 is a coefficient in the Jacobian of transition. On both 
sides of (1.15J there are inverse radia of correlation in the quartic and in Gaussian density measure, 
correspondingly 0 ). 

For the limiting point of the interval of densities we have: 


It?* — 77cl: 



(1.16) 


where 77 = ^ cr is a diameter of a particle (hard core), 77 * is the limiting value of density for the critical 
regime. 

In general, coefficients Dq and Go in J1.13} depend on density. In this work we take them as 77 = >} c , 
where T) c = 0 .13044 is the critical density determined from jl.llf by the condition A = 0. The Hamiltonian 
£(po) in 1.8 is obtained as a result of integration in EU over all variables pt, with the exception ofpc 0 


1 While investigating the characteristics of the integrand, the origin of the coordinates is transferred to the point A = 0 , 77* = 0 , 
and thus the change of p 0 may accept both positive and negative values (see figurejlj. 

2 “critical” means that within this regio n the re exists the RG symmetry. 

3 The variable po is introduced in E in 0 as a result of integration over rfrjv and summation over N of the expression 5 (po - 
Nl vT/V)), where N is the independent variable of the number of particles. In order to integrate over p0 in |l.8| we surely have to 
take the grand partition function. 
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The chemical potential p appears only in E po 
with the distribution of the number of particles. 
By definition 


PV 

InE =- and 

0 


. Thus, the variable p 0 is connected with the density and 


<3InE _ <31nEp 0 
d(fip) ~ dp* 


(1.17) 


Based on these two expressions we get the equation of state: 


P=p{ T,i 7 ), (1.18) 

connecting such variables as pressure, temperature, density, p, t, 77 . The curve of equations jl .17} is an 
envelope curve for pV/6 — InE, the latter depends on p, t, 77 , p. 

We work within the framework of the equilibrium thermodynamics. The number N is large. We use 
the presence of the thermodynamic limit: 


N 

N —* 00 , V —<• 00 ; — = const. 

1 / 


(1.19) 


In this connection, when we integrate over p 0 in jl . 8 } we use the steepest descent method. The integration 
in jl. 8 j was carried out and discussed in (HE). We give here only the final conclusions. 

In the integral /exp[/V£(po)]dpo in 1 1.8' we looked for the maximum of the subintegral function. 
The derivative 

d£(po) * o 

—- p + 2Dpo - 4Gpii = 0 (1.20) 

dp 0 

gives us the cubic equations: 

Po + Vpo + W - 0, where 

on condition of E{po) < 0: 


1 D 

V=-, 

2 G 


p 

W = -—, 
4 G 


2D- 12Gpg < 0 or Stpjf 3 *) 2 + V > 0. 


The equation jl.21| has three roots, depending on the sign of the discriminant Q: 


w 2 y 3 

Q =-+ — . 

4 27 


( 1 . 21 ) 


( 1 . 22 ) 


(1.23) 


Here, the first addend is positive, and the second is always negative. Therefore, the discriminant Q of 
equation jl.21) may be positive, equal to zero and negative. The case of Q = 0 separates the roots. At 
Q- 0 (WI2) 2 = - ( V /3) 3 , we ha ve three real ro ots, am ong them only the root p™ 3 * = 2 \J-W12 = (/p*/G 
meets the requirement 1 1.22 | 4 From equation jl.23 at Q - 0 we get 


p* =+G 



and thus 


max _ , /_ 

Po ~ ± \ on' 


12 D 
3 G 


(1.24) 


The roots ' 1.24 are the vertex of the rectangle on the plane p *, p 0 . Let us introduce the notation 


, 2 D 

a — Gb , b- \ -. 

w 3 G 


(1.25) 


Out of the rectangle we have the region Q > 0, inside Q < 0. In the case of Q > 0, equation jl.20f , 
or jl. 21 j, has only one root CQ 


n max _ _ „ 1 . 

Po - Pi - <7 [ u 


1/3 


P* = ±G\B1\, 


(1.26) 


4 Two other roots, P 2 and P 3 = - 5 \/p* IG do not obey 1.22 
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Figure 1. The curve of the generalized chemical potential p* as functions of p™ 11 *. On this curve, the 
integrand exp[JVi?(p)] in the region of (1.8} has maxima [solutions (1.26} for Q > 0 and (1.27| for Q < 0]. 
Here, a = G( V2DI3G) 3 ~ t 5/2v ; b = V2DI3G ~ t v,2 \ d = \/DI2G ~ t vI2 . The thick lines correspond to 
principal maxima of the function E{p* ,po), thin lines correspond to secondary maxima, and the very 
thin line is a line of minimum. The origin of the system’s coordinates is located in the point p* = 0, A = 0. 
The abscissa axis is along the line p* = 0. The ordinate axis is along the rectilinear diameter A = 0. 


where 



q-4 at 7 = 0 , 

q -1 at 7 = 1 , 


which corresponds to (1.24} at Q = 0. 

In the case of Q < 0, we have three real roots: 


* y + 2kn 

p = acosy, and Pi,2,3 = wcos- 


(1.27) 


n = 0,1,2. Only the roots 


<p w + 2 n 

pot = b cos — and p 02 = b cos —-— 


(1.28) 


meet the requiremets (1.22} . The principal maximum of exp[NE(po)] for the root P 02 takes place for 
| (the thick line in figure [lj, and for the root poi, when 0 cp ^ |. The secondary maxima of 

exp[NE(po)] exist for P 02 in the region 0 (p < | and for poi in the region | ^ cp ^ n. 

The curve for the principal maxima of the generalized chemical potential p* has a jump in po from 
—d to +d on the axis p* = 0 . 

A smooth transition of the solution p 1 into P 02 at the point (- b,-a ) and pi into p 01 at the point ( b , a) 
takes place. 

Herein below we shall consider an isothermal process of a quasistatic pressure on the gas system at 
r = const, r < 0 , [r| < t*. 


2. Monophasic states of liquid and gas 


The solution p^^ = pi in 1.26}, describes the values of densities of absolute maxima for the function 


exp[NE(p)] in integral (1.8} for E 
descent method, we write: 


Po . Considering A) as a very big number, according to the steepest- 
lnEp 0 = A7(l - A)p* + NE(pi). (2.1) 
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Eliminating p* by means of 1 1.17 



« in =„ =JV 

dp* 

(2.2) 

one gets (1 - A )N + Np i = N, and 

Pi = A. 

(2.3) 

Then, in correspondence with (1.20} 



p* - p* (A) = -2DA + 4GA 

3 and Pfp™ 351 ) = p*A + DA 2 -GA 4 . 

(2.4) 

Relationships (1.20)-(23) together with (T3)-(l3) give us the equation of state 



-P 0 V + P 1 V + Pp 0 V, 

(2.5) 

where 

PqV - 01nH o 

(2.6) 

— contribution from the reference system, 

P\V — 01nHi 

(2.7) 


is the contribution from the region ks= B and from the integrals over p k for 0 < k B. The latter includes 
integrals in the critical regime and in the inverse Gaussian regime, and the free energy of the RS up to the 
cumulant 9Jl 4 . 

Finally, 


P Po V = Q ln~ Po , (2.8) 

lnHp 0 = NS[A), (2.9) 

where in accordance with (2.1H2.4} 

<?(A) = p* (A) + DA 2 - GA 4 , (2.10) 

p* (A) = -2D A + 4GA 3 , (2.11) 

P po 1/ = QN (-2DA + 4GA 3 + DA 2 - GA 4 ). (2.12) 


The sum of (2.6} , (2.7} and (2.12 in (2.5} presents the gas-state isotherm when A ^ - fa, and the liquid-state 
isotherm when A 3= fa. 

Our main task is to consider events in the region |A| *£ b. In this region, the discriminant (1.23} of 
equation (1.21} is negative. Note that we restrict ourselves to t < 0, |t| <t* - 0,02 and 


|t?-t7cl < Be 



(2.13) 


By definition, the coordinates of the critical point are r = 0, Afpc) = 0, p c = 0.13044. 

The expression for A is given in l |l.ll} . The quantity A is also determined in (23). Here, for p \ we 
have its value (136), the quantities D and G being connected with the attractive forces. Thus, owing to 
expression s both the connection of the effects of the long-range and of the short-range (of the RS) 
interactions takes place. 

In the critical region, we can also use for A the following expansion: 


A (p) = A(/7 c) 



(■n-TJc)- 


(2.14) 


Here, the following expressions hold (dA/dp) J)=)? = According to the definition of the critical point 
A\ v=ric - 0, and for (2.14) we get (see figure (2) : 

A = — (77 - /yd (2.IS) 

be 
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Figure 2. The region of the linear approximation for A as a fun ction of the density and vice versa. Our dis¬ 
cussion concerns the critical region of densities [see equation 1.16 ] I 77 — 77 c! ^ 0.02 A = -('Df 2 f + 3 ®f 3 ^ 2 ), 
ri = [NIV) [ 71 a 3 16), rj c = 0.13044. 


and 


17 — T)c\ + T lc ■ 


(2.16) 


It follows that for gaseous state pi^-b and for liquid state pi 3= b in 2.IS ■, one has some restrictions 


on values of A and T). Besides, we should take into account the restrictions on density in {2.13} , as well as 
those connected with cumulants Ht 4 (2j[5j. 

For p* we have expressions {2.4} and {1.9} . The latter includes the chemical potential of the reference 
system, and has a contribution into the pressure of the whole system. 


3. On the reference system 

Let us turn back to formula ClD for the grand partition function H. The logarithm of H gives the 
whole equation of state accounting for external shell, thus accordingt to 1ft— — -^sph^sph ^Ssph> 
where according to O P S ph - p + SPq. In this section, we omit the terms connected with the existence 
of the external shell assuming InS = PVId. In accordance with dU 

InS = InSo + lnHi +lnH po . (3.1) 


Our aim is to separate in the right-hand side the terms belonging only to the RS. At first it is In —a - Pq V10, 
where Pq is the pressure of the RS. Then, from 2.9 , ' 2.10 and jl.9| for lnH po we may write: 


In 


H po = N^p-p 0 ) + S- |a(0)|(l - A)] + 7V(DA 2 - GA 4 ). 


(3.2) 


Thus, in the right-hand side we have two addends connected with RS, namely, N^-fip 0 ). They do not dis¬ 
appear, when excluding the attractive interaction, taking a = 0. Of course, in the right-hand side of {33} 
we should exclude some additional expressions, completing together the free energy Fo of the RS. We are 
talking about the factoi[^]lnQ -1 (-£, a 2 , a 4 ) that forms a part in the expression for Hi. Collecting together 
all terms remaining in the expression {3.1} at a - 0, we can write {3.1} in the form: 

PV P 0 V 

— - ^- = N{p-p 0 }-{F-F 0 ), (3.3) 

5 It is shown in O, Q(-f,« 2,04) = [/^exp|-fp; -y^- ^fp 4 ] dp;] B , where «2,«4 are the coefficients in the transition 
Jacobian into the phase space of collective variables. 
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where P 0 VI© = ln^. 0 . and 
F-F 0 


]V | -1 a (0) | (1 - A) + N (DA 2 - GA 4 ) ] + 


+ In Eg - /3 (F C r + -Figr) - N 



Nb 

N 


lnQ(f, a 2 ,« 4 ) 


(3.4) 


As a result, we may write Fq in the form 

^ = NS-N B \nQ(S,a 2 ,ci 4 ). (3.5) 

Herein below in our computations, in the critical region we shall consider the quantities entering 
{3.4) as the known quantities. We have presented here the expressions for Fq because we say that the 
pressure Pq, and the chemical potential po ought to be determined from the expression for Fq. The latter 
is obtained here using the quartic density measure while calculating the Jacobian. Thus, for the chemical 
potential of the RS we shall take the expression 


dF 0 

P ° = d^' 


(3.6) 


The extraction of the RS while calculating the equation of state according to {3.3) , will not be considered 
here, and will be the subject of another research. 


4. The field of the phase transition 


We consider the partition function HD in the case of negative values of the discriminant {1.23) of 
equation {1.20) . The part of the expression HD will be analysed in detail, which is connected with the 
chemical potential via E Po according to {1.17) . For E po , we have expressions {1.8) -{1.13). We shall remem¬ 
ber here the main form 


-Po = ex P 


[iV/i* (1 - A)] J exp[jV£(p)]dp. 


(4.1) 


In the region Q < 0, equation {1.20) for extremum of exp ATE(p) has the roots of {1.27) . 
We shall investigate the equation {1.17) for the envelope curve in the case of Q < 0. 
Let us consider more in detail the second of the two equations 


Taking into account {4.1), one has 


dlnS po 

dp* 


— N - 0. 


f Np exp(ATE(p))dp 

iV(l - A) + -= N, 

/ exp(jV£(p))dp 


(4.2) 


or 

OO 

/ pexp[iV£(p)]dp 


/ exp[NE{p)]dp 

—OO 

On the right we have A = A(p), given in {1.11) , where the cumulants of the RS are present. On the left- 
hand side of equation {4.3) there are the magnitudes connected with the long-range effects of the Van 
der Waals forces. This important equation connects together the long-range effects accumulated in the 
left-hand side, and the short-range effects accumulated in A. As we remember, long-range forces are 
described in the collective variables phase space and the short-range ones in the Cartesian-phase space. 
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Phase transition of the first order 


We devide the integral p exp[NE{p)]dp into two parts 


CXJ 

/ 


u 

/ 


p exp[ATE(p)]dp = / p exp[NE{p)]dp + / p exp[ATE(p)]dp. 


CXJ 

/' 


In the first one, along with figure[l] the maximum of the subintegral function exp[AT E{p)] takes place for 
p - P 02 - b cos 


<p+2n , 0 ^ tp 7 r. This is the domain of the gas state of the matter. In the second term, the 


maximum for exp [NE[p)] takes place for p — pot = hcos j — this is the liquid domain. 


In the whole integral |4.3) , the principal maximum for the gas-phase trips together with the secondary 
maximum for liquid, and vice-versa the principal maximum for the liquid phase trips together with the 
secondary maximum for the gas. 

Let us go to a more detailed consideration of expressions |4.3) . 

So, we carry out the integrals in l|4.3), and denote by 


J7 = 


= J p exp[NE[p)]dp, and by JE = J exp[NE(p)]dp. 


(4.4) 


Taking into account that the maxima of the subintegral function go along the lines p = po 2 and p = pot 
and that the curve P 02 , where —b =£ P 02 <-bl 2 , is connected with the gas-phase, and the curve poi, where 
b/2 sS poi b, is connected with the liquid phase, we shall write down integrals |4.4| as the sums: 

J* = J<g + A and = + (4.S) 

where 




U 

/ 


p exp[ATE(p)]dp, 


—OO 
0 


^g = 


/ 


exp[NE[p)]dp , 


OC 

I 


A- I pexp[NE(p)]dp, 


OC 

I 


JC 1= / expIA/iETpbdp. 


(4.6) 


(4.7) 


All of them are of the same type, and we shall calculate J? g precisely. We expand the function E(p ) at the 
maximum point p = P 02 - b cos <P+ ? 7T . Then, 


where 


Jg = exp[N E(p 02 )] (AiP 02 + A 2 ), 


(4.8) 


Ai = 


U 

/ 


exp 


-y|£(p 02 )l(p-p 02) 2 


dp = 


yjn 

— 


N .. 


-\E{p<a)\ 


- 1/2 


IP02 I 

/ 


exp 


N .. n 

--\E{p m )\x 2 


dx 


-^12 = 


u 

/ 


(p - p 02 ) exp 


- Y|£'(po2)l(p-po2) 2 


dpo 2 : 


777*:-- exp - - AT| £(p 02 ) [ p 02 ■ 

\NE(p 02 )\ 2 


For the second integral in Ai, we introduce the probability integral and ge{^| 


IP02 1 

/ 


exp 


N .. ? 

~—E(p 02 )x 2 


dx = (erfcz) 


\/f[£(P 02 )f 


6 While calculating Jhi, J 12 and JT we used expiv[.E(po2) - ||£(P02)Pn?l - exp[jV£( 0 )], £( 0 ) = 0 . 
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N - 11 __2 

Z = \ —\E{po 2 )\\po 2 \, erfcz = —=-e 
2 y 7T Z 


1+ E (-«' 

m=l 


1.3... (2m- 1) 


(2 z 2 ) m 

Substituting into {4.8} and gathering together all the terms after shortening in the limit of big N, we find 

\fn 


= P 02 e 


NE(pQ2 ). 


y/%mpta)\ 

The integral./] when calculated undergoes the same transformations and we obtain 

= p m z NE{pm) - 


(4.9) 


\fn 


\J f |£'(poi)l 


(4.10) 


The integrals J£ g and Jd\ are actually calculated while calculating - r /\ \. For Jd g and -7{\, we have: 


Jdg = exp[iV£’(po 2 )] 
= exp[AlB(poi)] 


x/tr 


^f[i?(po 2 )l 

s/n 

)/fl^(P0l)l 


(4.11) 


After all these calculations, we return to expressions {4.2) , {4.3} , that together with {4.1} describe the 
envelope-curve equation, or, in other words, to the equation of state in the region of Q < 0. Expres¬ 
sion {4.3} now has the form: 


7 ?g + ^i 

A = —-—, or explicitly A = pazWg + pot 


Jffg + 


(4.12) 


where 


nv 


c g exp[Al£(p 0 2)] 


JfTg 


mi 


c g exp[lV.E(po 2 )] + c 1 exp[A/'£(p 0 i)] ~£g + J£\ ’ 

_ qexp[Af£(p 0 i)] __ J7\ 

c g exp[NE(p 02 )] + Ciexp[ATE(poi)] J?T g + J£i ’ 

c g =([£(po 2 )ir 1/2 , q = [£(poi)l" 1/2 . 

Here, the functions m g and w\ play the role of the probabilities of the gas or of the liquid state, accordingly. 


As we see, 
Herein, 


For A, we have: 


W\ ■ 


1 

at 

<p - 0 , 

\ 0 

at 

<p= 

0 , 

1 

2 

at 

<*> = !- 

w g ~ 1 2 

at 

<p= 

71 

2 

0 

at 

cp = 7i; 

l 1 

at 

<p= 

71. 


P02 = £>COS 


Wg + W\ — 1. 

(p + 2n 


-b, cp-n, 
-d, (p — |, 
-f, ip - 0. 


CD 

Pot = & COS - = 


2 * <P = *. 

d, (p = |, 


b, (p-0. 
A = A s + A[, 


(4.13) 


(4.14) 


(4.15) 


(4.16) 
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Phase transition of the first order 



Pol 

Po 2 


Figure 3. Probabilities uq and Wg for single-phase systems, of liquid (uq) and of gas ( Wg ), respectively, as 
functions of the angle ip. Two bottom axes are to represent values pot (<p) and P 02 (<p) • 


) -b, cp-7i, 

- d/2, <p- f, 

o, <p = o. 


(4.17) 


{ 0 , cp-n, 
d/2, <p = |, 

b, (p- 0. 


(4.18) 


Here, b- \/|§ 


v/2. j_ ID _v/2 

’ “ ~ V 2 G _' 

In formulas l|4.13}, |4.17| and |4.18| we give important points for functions Wg, w\ and A g , Ap In fact, 
we should be driven by expressions (4.12) or even by more precise solutions of equation {4.3} . In figure[3] 
the graphics of Wg and w\ are presented as functions of the angle cp. 

Attention should be paid to the fact that the function A as a sum A = w g p 02 + w\p 0 i = A g + Aj turns 
into zero at A g = -d/2 and Ai = d/2. Infigure|4] the line {-d/2+ d/2) is interpreted as the line of the gas- 
liquid, first-order phase transition. The points A g | v=K /2 = -d/2 and Ai | v=jr /2 = d/2 do not lie on the curve 



Figure 4. Plot for A = A g + A] and A g = Wgpo 2 ; Ai = ifipoi; P 02 = b cos (p+ g ’' poi = b cos 0 S S i. 
The vertical fragment d/2 -4- -d/2 joins the point of the largest densities of gas with that of the smallest 
density of liquid in the two-phase system. 
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Figure 5. The quantities S, E and p* as functions of p max = b cos (fl+ ^ nn ,n = 1,2,3 within the region 
Q < 0 [see equations (1.23) , (1.27) and (1.28} ]. Thick lines represent the principal maxima for E{p). Thin 
lines represent the secondary maxima. For gaseous branch n- 1, the points -b,-d,-bl 2 correspond to 
values (p = 0,7iI2,ti; for liquid branch n = 0, the points b,d and bl2 correspond to the values ip = 0,jiI2,ji, 
respectively. 


A[(p)| v=7r / 2 . The sum of their values corresponds to the point A(</?) = 0 for <p = |. The expression A g = 
-dl 2 consists of two factors: P 02 = -d is the end of the principal maxima of the function exp[Ai£(po)l 
located upon the gaseous branch P 02 (</?), and w s = 1/2 that characterises the probability of arising of the 
value po 2 - For the point Aj = dl 2, it corresponds to w\ — 1/2 and the value poi = d, which starts the curve 
of the principal maxima of exp [AfE(p 0 )] located upon the liquid branch poi (</>)■ Note, that the density of 
either gaseous or liquid phase is independent of either poi or p 02 , but depends on A, according to (4. 16} - 
(4.18 '. A jump of the density at the phase transition is connected with the points poi = d and P 02 = -d, 
given in flgure[l]for p* = p*(po), but its magnitude depends on the values of products m g po 2 and w\pm 
in these points, in other words, on the probabilities w s (-d ) = \ and w\(-d) — | as well. Thus, we have 


Vg\ip=ni2 — p c[— d 12) + 77 c and m\<p=nl2 = Vc(dl2) + T] c 


(4.19) 


and the jump of density equals Tj c d. 

Also note that i?(p 0 ), which enter the coefficients c g and in 1 4.12 , are even functions of p 0 , -E(po) = 
2D- 12Gpg. Thus, E(b) = E(-b), E[d) - E(-d). In general, for the points P 02 and poi, which are symmet¬ 
ric relative to the rectilinear diameter A = 0, there is £(po 2 )/£(poi) = 1. 

In figure [ 5 ] there are presented the curves S = lnH po , E(p 0 ) and p*(po) as functions of p max = 
for n- 0,1, p m ax correspond to maximum values of the function E(p) in equation 1 4.1 for 


b cos 


( p+2nn 


InE 


pa- 


The very thin line represents the minima of the function E{p), which is the region of thermodynamic 
instability. The quantity exp[;VA(p) | characterises the probability measure for a state to exist. The points 
ff\ and / f g , which are the ends of f\ + ffi and / g 4- // g determine the limits for stable density fluctuations 
of liquid type arising in the gaseous phase and of gaseous type in a liquid. The points -bl 2 and bl 2 
determine the limits of secondary maxima for gaseous and liquid phases, respectively. 

The points -dl 2 and dl 2 are the points of fixed densities for gas in the gaseous phase and for liquid 
in the liquid stage in the two-phase system, which arise at the first-order phase transition. In figure [ 5 ] 
the closeness of the curves p* (p m ax) and &(p maK ) is obvious, which makes it possible, in what follows, to 
write <s?(A) ~ p*(A). 

Let us return to the initial expression for E po , given in (4.1| for Q < 0. Similarly to the case of Q > 0, 
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Phase transition of the first order 


in equations 12 . 8 1 -| 2.12 ' we write lnH po in the form: 


lnH po = NS [A), 


(4.20) 


keeping the notation of equations {2.8) - {2.12) . However, now expression {4.12) is taken for A. 

The quantity £{A), where A is a solution to equation {4.2) specified in {4.12) , is a curve that envelops 
the principal and the secondary maxima of the function E{p). 

It follows from figure [s] that the curves S{p max ) and p*(p m ax) stretch along each other. This is not 
surprising because <?(p ma x) consists of the sum p* (p m ax) and expressions Dp^ax - ^Pmax- the critical 
region p* ~ r 5/2v , and Dp 2 lax - Gp^ ~ t 3v . 

Therefore, in what follows, we shall write 


£ (A) = p* (A) + DA 2 - GA 4 , 
p* (A) = —2D A + 4GA 3 , 

£ (A) = p* (A) (4.21) 


and for the equation of state, taking into account l |4.20} and {4.21) : 

0 

Ppo = ^ ln -po> 

N 

P po =Q — p*{ A), -b^A^b. 


(4.22) 


The value of A is given in 1 4.12) . For p* we may take two equivalent expressions. In the first one, accord¬ 
ing to 11.27) 

p* = acostp, 0^(p^n, 

for the second one according to |1.20) 

P — _ 2Dp max + 4Gp^ lax , 

where, based on {1.27) , p max may take on two values: 

q) + 2n 

Pmax = P 02 = hcos---, 0 ^ <p *£ 7T, 


and this means the branch of the gaseous type, and 




Pmax — P01 — bcOS — , 0 


describes the branch of the liquid type. 


5. Properties of functions 8 (A) and A) in the region of Q < 0 


In the region of Q < 0 for 0 cp ^ n, there are two solutions to equations {1.20) , {1.21) : the solution 
of the gaseous type P 02 and that of the liquid type poi- Moreover, when po 2 (<p) for nl2 describes 

the principal maximum, in the same interval of values <p the solution poi(<p) describes the secondary 
maximum. Inverse situation is seen in the region where 0 ^ (p sS n/2: the principal maximum is due to 
poi, and the secondary one is due to po 2 - In general, this is governed by the probabilities w g and w\ given 
in {4.12). In order to take into account both maxima, in the expression for p* (p m ax) we should write: 


p* = p* (A(<p)) = —2DA + 4GA 3 


(5.1) 


where 


A = A g + Ai, 

Ag = w s po2, A\ = w\Poi ■ 
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Figure 6. The comparison of the curves pi* = acosip and pi* = -2DA[ip) + 4GA(y>) 3 as functions of the 
angle ip, 0 ip ji; N = 10 7 and N = 10 8 . As N grows, the attitudes of maximum and minimum remain 
unchanged, the widths go to zero. See also figures [3] and [4] for Wg and as well as for A g , Aj and A = 
Ag + Aj. 


Figure [6] presents the curves for the generalized chemical potential pi*{ip) as functions of A(<p): one 
curve is for pi* — acosip and the other is for pi* = -2DA{ip) + 4GA(<p) 3 . We are going to work with the 
latter one. They almost fully coincide except for the narrow interval of values ip near ip = nl2. It follows 
from figures[l]and[6] in the region of Q < 0, that in the process of isothermal compression of the gas along 
Po2(<f0 for ji 12 ^ ip sS n, the system via Ai enters the region of liquid-like values of pi* (A g + Ai). In other 
words, at compression of the gaseous phase in the region Q < 0, the density fluctuations of the liquid type 
emerge. Ultimately, at ip- nl 2, a droplet of liquid appears in the gaseous phase. 

What to do next depends on the type of the shell the droplet possesses when the two-phase system 
arises. The simplest case is when the shell of the droplet is presented as a geometric sphere confining the 
droplet. The shell does not contain particles. The question arises whether this is physically justified. In 
fact, the shell of the droplet has a volume, and is of the form of a spherical layer with a non-neglected 
thickness. All of them, i.e., two phases and the surface layer, should be present in the right-hand side 
of |5.1| . The left-hand side of {5.1} , which initially is of the form PVI& - lnH po , is written as a sum of 
terms dependent on the right-hand side of |5.1| , in the expression In H po . 

Returning to ED- for A (ip), in accordance with l |4.12} , we write 

A = A g + Ai, A g = Wgpo2, A] = w\poi. 

Then, 

S{ A) = <?(A g + Aj) = pi* (A g + Ai) + D(A g + A t ) 2 - G(A g + A[) 4 . (5.2) 

Let us write this as a sum 

8 (A g + Ai) = <fg(A g ,Aj) + <?i(A g ,Ap, 

pi* (A g + AD = pi* (A g , Ai) + nl (A g , Ai). (5.3) 

By doing so, we make no changes compared to either ED or ED- 

The curves for A g and Ai are presented in figure[3] Take the point A g = -dl2 on the curve A g , and the 
point Ai = d/2 on the curve Aj. Both correspond to equal probabilities w s = wi = 1/2 for Wg and mi in 
figure [3] Obviously, 

Ag{-dl2) + A\(d!2) = A[nl2) = 0. (5.4) 

In figures[l][s][6]the points A(?r/2) = 0 and pi* (A((p))A = o = 0 correspond to the transition from the gaseous 
part pi* (A(</?)) to the liquid part pi* (A(<p)). At this point, 

A = 0, pi* (A) = 0, <o(A) = 0. (5.5) 
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Phase transition of the first order 


Our aim is to show that the liquid-gas phase transition occurs along the line A g = -d/2 4^ Ai = d/2, and 
that condition (5.5) is obeyed starting at A g = -d/2 and finishing at Ai = d / 2 , and that this process can 


be described within the Gibbs statistics. Since the first-order phase transition at a constant temperature 
t occurs due to the external latent work of the pressure, we should demonstrate that at <p = n/2, and 
w\ = Wg = II2, the following is true 


Mg (<P)\<p=7il2 ■ 


: Mg (Ag.Ai) U g =—d/2 =0; 

A 1= d/2 


Ml (<M)[ip=7r/2 — Ml (A g ,Ai) a g =-d/2 — 0; 

A 1= d/2 

<§g((P)l(P=jr/2 — <Sg(Agi ^l) I A g =-d/2 — 0; 

A 1= rf/2 

\ip=nl2 — <ol(^ g i Ai) I Ag =-d/2 =0. (5.6) 

Ai=rf/2 

For this reason, we return to (5.2) and (53). Let us show that equations (53), (53) do take place along 
the line p * (A g , Ai) = 0 for A in the interval A g = - d /2 4- Ai = d / 2 . We write the expressions p* (A g + Ai) 
and £{ A g + Ai), selecting the mixed products A' g n A" and separating them for p g and pi and for S s and S\ 


p g (A g , Ai) = -2DA g + 4GA g + 12GA g A 1 2 , 
pf (A g , Ai) = -2DAi + 4GAf + 12GAiA g . 


(5.7) 


Obviously, the sum of these expressions is equal to the initial form of p* (A g + Ai) in ' 5.2 . Next, 


<?g(A g , Ai) = Pg(Ag,Ai) + DA 2 -GA^ + DAgAi 


-G 2 


i (A| Ai + A g Af) 


+3A 2 Af 


(g’lfAg.Ai) = pf (A g) Ai) + DA 2 - GAj + DA g Ai 


-G 2 




+ AgAf ] + 3A„A; 


2 a2 


g 1 


(5.8) 


The sum of these expressions equals <? g (A g , Ai) in 1 5.2 '. By substitution A g = w s p 02 and Ai = w\poi we 
make sure that on both ends of the line of phase transition 


A g = - d /2 - T - Ai = d / 2 , ( p - nl 2 , (5.9) 

(see figures[3]and[4) equalities (5.6) hold: 

Pg [- d / 2 , d / 2 ) - 0, pf (-d/2, d/2) = 0, 

(? g (-d/2,d/2) = 0, £\{- d / 2 , d / 2 ) - 0. (5.10) 

We conclude that the line - d /2 4- d/2 has the properties of the line of the first-order phase transition 
in the gas-liquid system, at arbitrarily large N. 

The expressions for generalized chemical potentials p * and pf as well as for £ g and £\ consist of both 
“pure” phase terms such as -D A g + 4GA|, and “mixed” products such as A g A 2 , A g Ai, A g Ai etc. On the 
lin e -d /2 4 - d/2, the “mixed” products are of the same order of magnitude as the “pure” ones. Note that 
in (5.8) , quantities p * , pf are proportional to r 5v/2 , and the expression DA 2 - GA 4 is proportional to t 3v . 

In figure[7j there are presented the curves for the mixed products A g A 2 , A 2 A], Their non-zero values 
are located near the limiting value ( p-nl 2 and characterise the width of transition layer between gas and 
liquid. The size of the region of values A g ”A[ ! depends on N. At N increasing, it narrows and the region 
of the surface shell joins the surface p * = 0, for the angle <p - n / 2 . Hence, the magnitude of the maxima 
keep their values unchanged. As a result, the effect of the mixed products on the requirements (5.5) , (5.6) 
remains unchanged. 
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Figure 7. The region of the surface layer described by the products A g A 2 , AgAj in the expressions for 
<?(A g , Ap and pi* (A g , Ap in 5.7 , 5.8 . The width of this region of values Ag'A” depends on N. At N 
increasing, this region narrows and joins the surface ip = n/ 2, pi* = 0. In addition, the magnitudes of 
maxima remain constant, and the effect of the mixed products on the requirement (5.5) , (5.6) does not 
change. 


It is seen in equation 1 5.7 that the mixed product 12GA g A 2 forces the gaseous branch pi * (A g , Ap into 
the point pi* = 0, A g = -dl 2, and the product 12GA g A 2 nullifies the liquid branch in the point pi* - 0, A = 
d/2, and we get the curve presented in figure[8] 

It remains to clarify the problem of the horizontal branch in the equation of state, the fact of originat¬ 
ing of the liquid phase in the gaseous one at isothermal compression of the gas, as well as the form of the 
left-hand side in equation (4.22) for the equation of state. 

This way, the existence of mixed products A" 1 A” in equations 1 5.7 1 , i S. 81 leads to the appearance 
of a shelf on the curve pi* ( A) = pi*{ A g ,Ai) + p*(A g ,Ai) in the equation of state. Conditions "" 
fundamental significance because condensation into the liquid starts at density of gas Tjg 
and terminates at density of liquid rj\ - 77 c(1 + dl 2). In other words, the gas isotherm ends at the point 


s 

= rhd 


have a 
-d/2), 



Figure 8. The curve of pi* (A) as a relult of combination of pi* = acosip and A = iv g po 2 (<p) + rr’ipoi (<p). The 
influence of the mixed products 12GA g Aj 2 on nullifying the branch pi* (Ag.Ap in the point A g = -d/2 and 
of 12GA g A 2 on nullifying the branch p* (A g , Aj) in the point Aj = d/2 is shown. 
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Phase transition of the first order 


P 02 = d, A g = -d/2 at density p g = p c (l - dll'], and the liquid-phase isotherm ends at the point p 0 i = 
d, Aj = dl2 at density r\\ = p c ( 1 + dl2). On the line p c ( 1 - d/2] + p c (l + d/21, the phase transition occurs 
from the gas into liquid. This phase transition from gas into liquid along p* (-d/2; d/2) -s- p* (-d/2;d/2) 
occurs due to the external latent work of the pressure. 

Therefore, according to s there are contributions into the energy density of the surface layer due 
to p* (A g , Ai) and pj* (A g , A[) 

[ju 1 *(Ag > Ai)] gurf =12GA g A?, (5.11) 

3 

p,* (-d/2, d/2) - p* (-d/2, d/2) = - -Dd. (5.12) 

s 2 


Since A g = w g po 2 , Ai = wipoi, and the probabilities u/ g and w\ depend on the number of particles N and 
the angle <p, as is seen in figure[8] At N increasing, the products of w g and w\ degenerate into 5-functions 
placed along (p-nl 2. This means that in the limit of large N, the two-phase liquid-gas system degenerates 
into the system gas — a droplet of liquid divided by the surface without particles. 

The energy of the phase layer transforms into the energy associated with the creation of a geometric 
surface of division. The two-phase system emerges: gas and a droplet of liquid of radius R<i r . If for the 
gaseous phase P g is the pressure, V g is the volume, P g V g is the energy, then for the droplet of liquid the 
volume is ^R^ r , the pressure inside it Pi = [P g + 2a/i?dr). where a is the surface tension coefficient, the 
energy P\V\ — ( P g + 2n//?dr) / 3) and plus the surface tension energy -aSdr = -a4nR^ r . 

Hence, the formation of the droplet with a geometric interface is concerned with an additional energy 
of the surface 


2a 4 nRl r 
/'i;r 3 


- aAnR 2 


4ti R 


dr 


Rdr 


(5.13) 


where the surface-energy density is -a/i?dr- The magnitude of the surface energy is connected with 


quantity (fsTTl). Also, the change of the volume for the gaseous phase takes place, from V\ to 
(4nR? I3){n\l n g ). 


Va - V - 


6. The equality of the chemical potentials of gas and liquid on the 
condensation-boiling line /i*(A g , Ai) 

For p*, by definition, we have equation |1.9| : 

p*=/3(p-po) + <-|a(0)|(l-A). (6.1) 

In the approximation accepted in this research, the chemical potential of the reference system is 

|Sp o = £. (6.2) 

Then, for the chemical potential p we get the following expression: 


p = -0|a(O)|(l - A), 
A = w g p 02 + wipoi. 


On the line of phase transition the conditions 1 5.5 i, 1 5.6 1 , p* (A g , Ai) = 0, hold together with 


N g + N\ = N, dA/g = -dATi, 


(6.3) 


(6.4) 


where N g and N\ are the number of particles in the gas phase and in the liquid phase, respectively, p g , pi 
are the chemical potentials of the liquid and of the gas, respectively. According to j5.6} 

j8pg = |o(0)| g (l-A g ), 

)Spi = |n(0)[i(l - Ai). (6.5) 
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According to definition (1.12) 


N 1 - 

|a(0)| =V0 I ® (O)I - 

|a(0)| g = J7 g ——~ ^|®(0)|, 

6 b no A 0 

l«(0)li =171—3 ^l®(0)l. 
no A 0 

We consider the beginning of condensation (boiling) 

A =A g + Ai, 


then 

p g — gas density, 

T)\ — liquid density in a droplet. 


(6.6) 


Ag-“P 02 > P 02 - -d, Ag -~-d, 

Ai = -poi, Pot = d, Ai = -d. 


(6.7) 


The densities rig and tj\ are expressed via the value ? 7 C in the critical point and via A based on the rela¬ 
tionship j2.15| : 

77 = 77 c A + ?7 c, J7 g = 77c[l--d], m = 7 7c(l+ -d). (6.8) 

Then, for (5.4) we obtain 


Pg = l®(0 )|—3 77c 1 


71(7° 


M 1 = I $( 0 ) 1 —o* 7 c 1 - 




=(> 4 )- 

(■-?)■ 


(6.9) 


Thus, 

Pg-p\. (6.10) 

It was proved that on the line p* = 0 there exists a two-phase system of a gas phase and of a liquid phase: 
the maternal sphere of the gas and the droplet of liquid inside it and the chemical potentials of both 
phases are the same. Let us imagine that conditions (6.7) and ( 6 . 8 ) necessary for equalities (6.9) and (6.10) , 
may be also valid for other quantities A g and Aj, not only for that of (6.7) , which are located symmetrically 
to the point of the intersection of the rectlinear diameter with the line p* = 0 that is A g = -A], Such points 
are located at the ends of the intervals which pass through the origin 77 * = 0 and touch the lines p 02 and 
pot that are the curves of the maxima for gaseous and liquid phases, respectively, see figure [T] For each 
pair of points, denoted by po 2 and poi, the equalities f?(po 2 ) = E{p 01 ) take place, since P 02 = -poi. and 
£(po) is an even function. Thus, for each pair of end-point values there is w{poz) = m(p 01 ). Hence, for 
each pair, the following is true: w{p m )P 02 = ~w(.poi)poi and A = m(p 02 )po 2 + w(poi)poi = 0, p* = 0, 
S{ 0) = 0, and also p g (A g ) = pi(Ai). However, conditions 1 5.6 1 are not obeyed if A g = -d/2, Ai = d/2 are 
substituted by any of A g = -K\. 

Therefore, neither approaching nor moving away from the densities of gas and liquid compared 
to (6.7) occur. The conclusion is that densities for gaseous and liquid phases do not change at the first- 
order phase transtition. The density of gas remains fixed 


Vs 


■ p c Ag + ??c = 77c(l - d/2j. 


The density of liquid remains fixed 

771 = 77cAi + 77 c = 77c(l + <^/2j. 
The chemical potentials of both phases are the same 


( 6 . 11 ) 


( 6 . 12 ) 


p g = pi. 

This remains true as long as the phase transition takes place. As a result, under the action of latent work 
of the pressure, all the gas condenses into liquid. A jump-like change of the gas density takes place 

Pi-Pg-flcd. (6.13) 
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7. The process of condensation of gas into liquid (the boiling of liquid) 


This section begins with (4.22) 

P V 6 

-= NS {A), P = ©ng( A), 0 = 0—j. (7.1) 

0 TUT 6 

Let us follow the process of isothermal quasistatic compression in the gaseous phase (see section]!). The 
process takes place along the line A = w g po 2 + uipoi- Starting at A = A g = w g p~o 2 : 


Po2 = ~b, w g = l, A --b, w\ — 0, (7.2) 

where (see figures]!]]!) there is a homogeneous gaseous system (see figures[l]and[5). 

Next, after the quantities p 02 and A increase, we get p 02 > -b, A > -b, and the probabilities are as 
follows: 

w g < 1, mi > 0, Wg + w\ - 1. 

In a gaseous system, there emerge density fluctuations of a liquid type. During the quasistatic process for 
-b « P 02 ffi fluctuations are unstable, because, according to figure [s] under the line f f\ - f\ E(A) < 0 
and exp [AH: (A)] —• 0. 

Being further compressed, while crossing the line (see figure ]s) in a gaseous system, there 

emerge stable density fluctuations of a liquid type because E( A) >0 therein. The argument A is deter¬ 
mined by equation (4.12). Here, both terms in A become of non-zero value. However, the term w g po 2 in 
A is a principal one. As earlier, we have only gas, and there is a single-phase system (see figure]!). 

When the external pressure increases, the pressure in the gaseous system also increases. The den¬ 
sity of a gaseous system grows, in accordance with (2.15) . The gaseous system evolves along the same 
isotherm p* (po 2 ). Now we get to the limiting densities p g = p c ( 1 + Wgpo 2 ), rp — ri c [l + m;poi), where 

P 02 — —d, poi-d, E{d)-E{-d), Wg-w\-l/2. (7.3) 

Based on (4.12) , (5.6) , (5.8) - (5H0) , it follows 


A = A g + Ai = 0, p* — p g + pj* = 0, § — Sg + S i = 0. 


(7.4) 


This state describes the end of the curve for the principal maxima for a gaseous phase and the beginning 
of the principal maxima for a liquid phase. Here, we arrive at (5.10) and m 

When work is done on the gaseous system, within the gas system there arise processes that evolve 
along the curve Ai = wipoi, although under the conditions of a secondary maximum, as is seen in figures]!] 
and[s] The system does an internal work to create stable density fluctuations of liquid type, which in the 
case of p 02 = —d, A g = -dl 2, poi = d, A\ - dl2 causes the creation of a droplet. 

In order to determine the properties of a droplet, it should have exact boundaries. Equations (5.7) 


and 1 5.8' for p *, p*, pj\ <?, <? g , g\ are “bulk” formulas. In these formulas, the interface is a phase layer 
with a number of particles. Our further calculations are performed for the case of very large N. The 
model of a two-phase region is considered as gas and a droplet of liquid arising in gas, the interface 
is a geometric surface without particles. The mixed products A g n A", entering p* (A g , A]) and <?(A g ,Ai), 
presented in and (5.8) are replaced with the energy of additional pressure in the droplet of liquid 
(4n:Ji^ r /3)(2a/if t ir) an d the surface energy -a47iR^ r . 

The two-phase system exists only along the line p* = 0. On this line, the chemical potentials of gaseous 
and liquid phases are the same. Equality (5.11) is invariant of the process of condensation (i.e., boiling). 
This is true for the case when the thermodynamic potential <J> Po from the beginning to the end of the 
process of condensation (i.e., boiling) remains constant 


<h po = NgPg + Nipi = Npg = Npi, (7.5) 

since N g +N\ = N remains constant and is equal to the initial number of particles, and dlV g = - dN\, which, 
in fact, holds. The amount of gas decreases, while the amount of liquid increases till all the gas transforms 
into liquid. 

Liquid phase originates in the form of a droplet confined by the interface S of some surface energy. 
The sharp bend of the isotherm in figure[8]is related to the instant of the creation of a droplet, and the 
sum of the moduli of the mixed products (5.11) is equal to the energy of the droplet. 
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7.1. The energy of the surface tension 


According to ' 4.12 <, |6.5|-|6.7|, along the line /i* = 0 we have, A = A g + Ai, A g = w g (-d) = -dl 2, Ai = 


injd = dl2 , u> g - w\- 1/2. 

For the basic relationship j7.1), the left-hand side takes the form 


PgVg + Pi Vi — a Si, 


(7.6) 


where P g , Pi are the pressures in the gas and in the liquid phases, V g and V\ are the volumes of the gas 
and of the liquid phases, Si is the surface of the droplet of liquid, a is the surface-tension coefficient. In 
the right-hand side of (Z3 we obtain: 


N g S g (\ g ) + - Afg/i* (A g ) + N^* (Ai). 


(7.7) 


According to the conclusions drawn in ' 6.5 and |6.6| , the quantities /i g , p* take on fixed values. In accor¬ 
dance with d6.5M6.7l: 


<gg(A g ) = n* (A g ) = -2DA g + 4GA 3 = 3Dd/4, 

<g\(Ai) = /i*(Ai) = -2DAi + 4GAf = -3Ddl4, (7.8) 


and the process of the phase transition is reduced to redistribution of N g and N\. The width of the inter¬ 
face between gas and liqud droplet is considered to be zero. 

Liquid phase is a droplet of radius R that originates in the gaseous phase. The droplet has strict bound¬ 
aries. Combining j7.5) - j7.8| , for iz3 on the line fi* - 0 we obtain: 

P g Vg + Pi Vi - aSi = &{N g S g + NiSi). (7.9) 


The emergence of a liquid droplet, restricted by a shell, in a gaseous phase is the essence of the first order 
phase transition. The droplet is of a spherical form. The surface tension energy is minimal. For j7.9) , we 
get: 


4n c, pi 

PgV g = PgV-P g -f4A 
3 'Is 

I 2a\4n~ 

m = \p.+^\—R 


h’dr 


dr ’ 


(7.10) 


where V is the initial volume, f?dr is the radius of the first droplet, rj i, p g are the densities of gas and liquid, 
respectively. For the energy of the droplet surface, we have: 

P 3 a 4 ttP 3 

-aS- -a47iRfc = — —--. (7.11) 

i?dr 3 

In the right-hand side of j7.9), it is written 


Q[N g S g + NiSH = Q[NS g + NiiSi - S g )]. 
Collecting {7.8M7.10) in equation j7.7) , we obtain: 


(7.12) 


P g V g + P l Vi-aS = P g V- 


/'dr 


4 ^dr n 4 ^dr 
^ g 


—-l 
Vs 


V — 


4nRl r 


— — 1 
Vs 


4nRl 


^dr 


= N&s g +msi - <? g ). 


(7.13) 


The change of gas volume by ^P 3 r (pi/p g - l), that takes place in the product P g [V - ^-R^ (pi/p g - l) ] 
in j7.13| corresponds to those changes in the quantity d? g (A g ) which were observed at the transition from 
^g(Ag) at Ag = -b to <og(A g ) at A g = w g {-d ) = -dl 2. We compare 


V- 


4nR i 


Vs 


= NQS g . 


(7.14) 
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Hence, for the specif] resurface tension energy of the droplet of radius R, we have 

~ — Yl\[§\ — Sg)©. 

K 


(7.15) 


Let us replace 8\ and Sg with 1 7.8and n\ with n\ = N\!(4nR^ r /3) - T]\6lna 3 . In the right-hand side of 


47.15) , we obtain: 


90 

Qn\{ 8 i- 8 g) - —771 — 5 Dd, 

° 7T/T-J 


(7.16) 


which conforms to (5.10) . Therefore, based on 47.15) , there is a contribution to specific pressure P t otai 
from the surface tension energy of the liquid droplet of radius R: 


a 3 6 3 - 

ftotai = -- = —-O— 3 pi Dd=--®rnDd, 
R 2 tut 3 2 


(7.17) 


where 0 = 0 61 tut 3 . 

The quantity P t otai is proportional to r 5/2v . It governs the rate of decreasing the surface tension energy 
at approximatiion to the critical point r = 0. 


7.2. Latent work of condensation (latent heat of boiling) 

Condensation begins at the moment of the creation of the droplet with the surface energy (7.17) . As 
was mentioned above, condensation of a gas starts from the state 

NS {A) = JVp*(A g ), P02 — —d, A g = -^rf, p g = p c -^rfpc 

and terminates in a liquid system in the state NS[A) - Af/ij* (Ai), poi - d, Ai = \d and pi = p c + \dri c . At 
the moment when condensation starts, the pressure is equal to 


Pg = rig© 8 g{Ag) - rig©/!* (A g ). 


(7.18) 


Replacing A g with A g = w g {-d) - - \d and 


Then, 


1 , 3 

u*(A g ) = dD--d 3 G- -Dd. 
P g 2 4 


Pg = rig©-Dd. 


(7.19) 


(7.20) 


The initial volume of gas is a sphere of radius f?o, of volume Vo = 4ttRqI3. The final volume of gas 
after transferring into a liquid during condensation at a fixed pressure is the volume of liquid V\ in the 
form of a sphere of radius R\: 

n 

All N particles of a gaseous sphere Vo passed into a sphere of liquid V\. Thus, 

1 1 Vi 

Vo = V], Vi= —Vb, 

v g m vg 


the total change of the volume is 


(VS - Vb) = 1 - 1 1 Vb. 


The work performed is the latent work (i.e., heat) of the process of condensation: 


(7.21) 


-P g M-V 0 ) = -P g -i-1 \—R?. 


v\ 


4 n. 


'By a volume unity. 
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Having substituted 1 7.20the specific latent work of condensation i^] 


L = - |l - y j n g @^Dd = -{v g - i^©Ddn 2 g . (7.22) 

Condensation at isothermal compressibility is a dynamical process. Equation (7.9) concerns the states 
of a system in which two phases have already been created: the initial gaseous phase and a new liquid 
phase in the form of a droplet. The conditions A = 0, p* = 0, ^(A) = 0 require an external pressure to 
be applied in order to do the latent work of condensation. As a result, the gaseous phase disappears (or 
remains as a nascent bubble of gas), and we write: 


P\V\ - L(V g - Vi). (7.23) 

[For a nascent bubble of gas, if any, P g y g nascent - aS g nascent is energy, and we can write: Pi Vi + P g u g nascent - 
aS nascent = i( y g _ ^nascent _ ^ 

For the two-phase system along the line p* = 0, S'{A) - 0 for external pressure we write down: 

^y=o = i(-^-l). (7.24) 

Using L, we get 

Pn*=o = \l~^n g ^©Dd. (7.25) 

On the curve of the equation of state, this expression represents a horizontal fragment of the transition 
from a gaseous phase into a liquid phase under the isothermal compression of a gaseous phase. 


8. The equation of state 


Now we shall summarize the results of our calculations in the critical region of r < 0, |t| < t* and 
I 77 — 77d < Ip* - Pel. where t* and |p* - p c | are the boundaries of the critical region presented in expres¬ 
sions (TlS) , (TTg) and (03) . 

The condensation of gas at isothermal compression is a process analogous to the boiling at an isother¬ 
mal decrease of pressure of a liquid. This is a quasiequilibrial directed process. Equation (7.9) concerns 
the states of a system in which two phases have already been formed: the initial gaseous phase and a 
new liquid one, in the form of a droplet. (The droplet has got a shell that separates liquid from gas, and 
its creation is the essence of the first-order phase transition.) Conditions A = 0, p* = 0, S’(A) = 0 require 
an external pressure to be applied in order to do the latent work of condensation (7.22) . As a result, in 
(7.9) the gaseous phase disappears (or remains as a nascent bubble of gas). 

The conditions existing on the line of condensation p* = 0, S( A) = 0, A g = -dl 2 + Ai = dl 2, p g = 
p c (l - dl 2) 4- p c (l + dl 2) = p/ are determined by (6.9) , i.e., by the equality of the chemical potentials 
of different phases. The thermodynamic potential <P po of a system on the line of condensation remains 
constant 

<P po — A/gPg + A/ipi — A/pg — fVpi, (8.1) 


according to (6.9 


O po = 7V|O(0)|—^-ped-rf 2 /4). 


( 8 . 2 ) 


On the line of the first-order phase transition, equations (7.9) hold. To obtain these equations, the mixed 
products A™A" were nullified in the exact relationships |5.7|-' 5.8 , in particular those in 1 5.11 1 . Instead 
of them, the surface tension energy of a droplet was introduced as well as an additional internal pressure 
of a droplet. This allowed us to describe the magnitude of energy associated with the creation of the liquid 
droplet. The droplet has got a strictly defined interface. We have described the quantity of the latent work 
of the pressure. 


8 Dimension of L coincides with the dimension of pressure. 
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Figure 9. The graph of <?(Ag) = p* (A g ) and <f(Ai) = p* (Ap; circles denote densities of the gaseous phase 
r/g = 77 c (1 - d/2) and of the liquid phase pj* = 77 c (1 + d/2), which remains constant during the liquid-gas 

phase transition. Both solutions are in the region of instability po 3 = bcos ( se e figure s[l] and [sj. 

The jump of the density at phase transition is 771 - p g = rj c d. The slope of the gaseous branch A g while 
approaching the point Ag = -d/2 is equal to 0. The slope of the liquid branch while turning away from 
the point Ai = d /2 is equal to 0 . 


The task of our further research will be to clarify the problem of the thermodynamic limit and to com¬ 
pare the theory with experimental results. We will need a total expression for InH, not only InH po . These 
urgent problems are easier to consider farther from the critical point, since there are no requirements 
following from the renormalization-group symmetry near T - T c . 

Thus, the general expression for the equation of state takes the form 


Ppo — 0 |t?gftg (Ag, Ai)^ g ^-rf/2,p*<-/!/2 

+ [“shelf”, the region of constant pressure]_/ ! / 2 cp*</ J /2,-d/2cA«d/2 

+ t?i/h (Ag, Ai)^ 9 d,p*>fi/2} > (8.3) 


where 


ft* (A g , Ai) = -2DAi + 4GAf + 12GAjAg, 
ft* (A g , Ai) = -2DA g + 4GAg + 12GA g Af, 


h —>• 0, 


h —>• 0. 


On the “shelf”, the external pressure is equal to the specific work of condensation according to (7(9), 
P p * = (l - vil Vg)ri g QjDd, and the thermodynamic potential remains constant, <I> Po = Np g {-dl2,dl2) = 
JVfii (-d/2, d/2). 

The phase transition is a dynamical process. It takes place without any breaks as a quasistatic isother¬ 
mal compression of a gaseous phase. There are five stages of the process. The first one is a gaseous phase 
with density 77 = T) c [-b) + T) c . Here, w s (q>) = Wg(n) — 1, u>\ = 0, A = A g = -b, b = \/2D/(3G), Ai = 0. 
When moving along A = A g = ifg((p)po 2 (<f>), there are density fluctuations of a liquid type in the gaseous 
state along the secondary maximum Ai = wipoi- The region of a sharp bend of the isotherm ft* (A g , A]) 
is reached due to the terms in jS.llf (see figure [9). In a gaseous phase, a liquid droplet appears with 
an additional, in comparison with the gaseous system, energy density J7.15| - f7(T8) . A liquid droplet that 
emerges means the end of the first stage and the beginning of the second stage, being the main stage of 
the first-order phase transition. 

The thickness of interface, accordint to figure [7] is equal to half-width of the maximum for the ex¬ 
pressions 12A g A 2 , 12AiA|, and 2A g Ai. In general, the interface includes all the mixed products 2A g "A 1 ,z , 


via ft* (Ag, Af) and <f(A g ,Ai). All terms of the mixed products distribute accordingly in equations ' 5.7 
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and 1 5.8' between £„ and 8\ in such a way that in the fixed limiting points of a phase transition A g = -dl2 


and Ai = dl2, the fundamental equalities |5.6| l hold. 

The third stage is considered to be the result of the action of the external latent work of pressure 
on the two-phase system, presented in j7.23) - fL24| i. Pressure and temperature do not change, accoridng 
to j6.5} - |6.6) . Gas transfers into liquid according to The external latent work of pressure has been 
done. The density jump Ap = rji - rj g - T) c d takes place. Everything ends with the fourth stage, which is a 
liquid with a nascent bubble of gas. Finally, the fifth stage is a liquid of density growing from the value 
771 = Ti c dl2 + T) c to 771 = Tjcb + r/c. Thus, the events in the region Q < 0 finish. The region of liquid densities 
begins, where the discriminant Q > 0. This is a homogeneous liquid phase. 

We started with an initial sphere of gas in a boundless system. By isothermal compression, we reached 
a point of the first-order phase transition. Gas transferred into liquid, condensed into a droplet. The 
quantities that describe the process were determined. We worked within the framework of the grand 
canonical distribution for the Gibbs statistics. The successful solution of the problem is ensured by the 
choice of a related phase space of collective variables, among which there is a variable po- This variable 
is responsible for a description of the phase transition and the density jump. Calculations are carried 
out using the quartic basic measure density. The zero-approximation performed in this work provides 
good results on both the critical indexes and the description of the mechanism of density jump at the 
condensation of gas into liquid. 

We have gone through all the events of the first-order phase transition below the critical point. The 
curve for the equation of state has a form presented in figure [9] A peculiar thing is the density jump 
between the points p g - T] c ( 1 - dl2) and 771 = tj c (1 + d/2), where d - VD/2G, where D and G are the 
coefficients at the second and fourth power of p 0 , respectively, in the initial Hamiltonian E(p 0 ). 

Figure |T] resembles the Van der Waals isotherm. The theory put forward by a great physicist fTTi , 
was the basis for numerous researches devoted to the gas-liquid phase transition, including the present 
research. 

At the same time, this work is not void of numerous evident drawbacks. Among them we note: 


• we worked in a very narrow regions of temperatures and densities, near the critical point; 

• when calculating the Jacobian of transition to the phase space of collective variables, the statistical 
weight was introduced, the distribution over short-range interaction between hard spheres. The 
cumulants of the Jacobian, SRi, W( 2 , 'JJi-. and 9Ji 4 , were taken at k = 0. This was connected with the 
fact that they have wide “shelves” at k = 0. The cumulants were determined via compressibility of 
the reference system. We neglected the terms proportional to k 2 , which were proportional to the 
binary correlation function for the system of hard spheres; 

• we neglected the effect of integration over p^ for k + 0 on the form of the isotherm of phase transi¬ 
tion. 

• we restricted the expression for the chemical potential of the reference system to /3po = 

We postulated that the dependence £{p m ax) = (A), being accurate in the region of Q > 0, remains 

also accurate for Q < 0, if the sum A = A g + Ai is assumed as A. 

All the above demands a more detailed consideration. The mutually-coordinated short-range and 
long-range interaction being taken into account will be one of the main tasks in further research. 
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0a3OBi/m nepexifl 1-ro pofly b 06 /iarri kpmthhhoi tohki/i 
ra3-piflMHa 

I.P. K)XHOBCbKI/m 

iHCTMTyT tf>i3HKM KOHfleHCOBaHnx ci/icreM HAH YKpaiHW, By/i. I. CBeHqiqbKoro, 1, 79011 /IbBiB, yitpama 

Po3r/iflflaeTbai noBefliHKa ci/icreMi/i B3aeMOfliroHi/ix HacTi/iHOK b o 6 /iacri TeMnepaTyp HWXHe KpwTHHHOi tohki/i 
T sS T c . 3aBepLuyeTbcti po3paxyHOK Be/utKOi CTaTMCTHHHOi cyMW, notaTi/iii y nonepe^Hix po 6 oTax y Merafli ko- 
zieKTUBHUx 3MiHHHx. 3a 6a30By rycTHHy Mipu 6 epeTbcn neTBipHi/iii (a He TaycoBHii) po3nofli/i. OnncaHi noflii 
noBh3aHi 3 (J>a30BMM nepexoflOM 1-ro po/ty, mo Bifl 6 yBaK)Tbcn b pe3y/ibTaii i30TepMi4Horo KBa3iCTa™HHOro 
CTHCHeHHti ra3y. Bnfli/ieHa /liHin /.i * ( 77 ) = 0, Ha rk\v\ y ra30Biw <f>a3i ni/; flieio Ti/icxy bhhhkbc Kpan/in pifli/iHH. Mae 
Micqe piBHiCTb xiMiHHHX noTeHqia/iiB ra30BOi i piflKoi' (y Kpan/ii) c)>a3; 3Ha(ifleHO Be/iHHHHy noBepxHeBoi eHep- 
rii xpan/ii, po3paxoBaHO cxpHTy po 6 oiy KOHfleHcaqii, Bi/i3HaHeHO cxanoK rycTHHH, HamtcaHO piBHnHHn CTaHy. 
Po 6 oTa CTaHOBMTb neBHy KmqeBy CTafliio floc/iiflxeHb b 06 /iacri TeMnepaTyp i rycTHH, m° BK/uonae y co 6 i Kpn- 
TMHHy TOHKy piflHHa-ra3. 

K/itoHOBi c/ioBa: KpuTum-ia ro^Ka, <pa30BMi/i nepexip, nepinoro pop,y, Beni/ixa CTari/iCTi/iHHa cyMa, KOJieKTi/iBHi 
3MIHHI 
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